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We present a Monte Carlo study of a lattice gas driven out of equilibrium by a local hopping bias. 
Sites can be empty or occupied by one of two types of particles, which are distinguished by their 
response to the hopping bias. All particles interact via excluded volume and a nearest-neighbor 
attractive force. The main result is a phase diagram with three phases: a homogeneous phase, 
and two distinct ordered phases. Continuous boundaries separate the homogeneous phase from the 
ordered phases, and a first-order line separates the two ordered phases. The three lines merge in a 
nonequilibrium bicritical point. 
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I. INTRODUCTION 

The statistical treatment of systems driven far from 
equilibrium presents exciting theoretical challenges [1, 2]. 
Lacking the unified understanding afforded equilibrium 
phenomena by the work of Boltzmann and Gibbs, we are 
exploring unknown territory without recourse to an es- 
tablished theory. Our physical intuition, developed in the 
context of equilibrium systems, can be misleading when 
faced with nonequilibrium problems. We therefore turn 
our attention to computational 'experiments' in which a 
manifestly nonequilibrium state can be established and 
studied, seeking to identify key features shared by many 
nonequilibrium systems. While real physical systems are 
unquestionably important, complications and subtle de- 
tails may obscure such common features. While moti- 
vated by realistic problems, microscopic rules and bound- 
ary conditions are chosen simple enough to facilitate a 
comprehensive computational study of the full parame- 
ter space. Many of the questions relevant in equilibrium 
remain interesting, especially as concern the nature of 
phase transitions and the principle of universality. The 
richness of nonequilibrium phenomena is often surprising, 
as the relaxation of the detailed balance constraint allows 
a variety of unexpected possibilities: in contrast to equi- 
librium, the dynamics now affects the stationary (long- 
time) properties of the system. Particularly dramatic 
effects have been observed in models where the violation 
of detailed balance is combined with spatial anisotropics 
and dynamic conservation laws [1]. There, effective long- 
range interactions can be induced even if the microscopic 
rules are perfectly local in space and time [3]. 

In this paper, we consider a model from this class, 
namely, a lattice gas of two species of particles and 
holes on a fully periodic lattice in two spatial dimen- 
sions. To drive the system out of equilibrium, we bias the 
hopping rates of the two species in opposite directions, 
reminiscent of an 'electric' field, E, acting on opposite 
'charges' (though we stress that there is no Coulomb in- 



teraction). A nonzero charge current signals the nonequi- 
librium steady state. The two species interact through 
an excluded volume constraint and nearest-neighbor at- 
tractions. We choose the interactions carefully, in order 
to unify three important models which appear as limit- 
ing cases of our more general theory. First, by letting all 
particles attract each other, irrespective of their identity 
(charge) , the non-driven limit corresponds to the familiar 
Ising lattice gas [4]. This well-known equilibrium model 
will serve as an anchor for our studies of driven systems. 
Turning the bias on but removing all members of one 
species, we recover the driven Ising lattice gas introduced 
by Katz et al [5] (the KLS model). The third limit, ob- 
tained by letting the interaction strength vanish, corre- 
sponds to a non-interacting two-species model first pro- 
posed by Schmittmann et al [6] (the SHZ model). While 
the KLS model phase-separates, via a continuous transi- 
tion, into high- and low-density strips aligned with the 
field [5], the SHZ model orders into transverse, charge- 
and mass-segregated, strips [6], similar to jamming in- 
stabilities in traffic models [7]. If a charge imbalance is 
imposed, these strips drift [8]. Our study will allow us 
to bridge the gap between these very different scenarios. 
To set the scene, we briefly survey some of the relevant 
previous work. 

As a single parameter modification of the Ising lattice 
gas, the KLS model is a minimal model for the study 
of nonequilibrium steady states (NESS). Particle hops 
along one direction (parallel to the z-axis) occur at the 
normal equilibrium rate, as if in contact with a heat bath 
at temperature T. Particle hops in the other direction 
(parallel to the y-axis) are enhanced (suppressed) in the 
positive (negative) direction by coupling to an external 
field, E. With periodic boundaries in the y-direction, a 
nonzero current is maintained, and the system settles into 
a NESS. At half-filling, there remains a continuous tran- 
sition, though with T c (E) increasing monotonically with 
E and saturating at 1AUT C (E = 0) [9-11]. The transi- 
tion falls into a novel universality class with exponents 
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distinct from the Ising ones [12]. The critical behavior 
is strongly anisotropic, with distinct sets of exponents 
characterizing fluctuations perpendicular and parallel to 
E. There has been some discussion regarding the nature 
of these fluctuations, with some authors disputing the 
original claim that the correct mesoscopic description is 
gaussian in the perpendicular direction [13]. Though the 
anisotropy makes numerical investigations of the criti- 
cal behavior quite subtle and computationally intensive, 
recent high precision Monte Carlo studies compare the 
two mesoscopic descriptions. The results are in complete 
agreement with the predictions of the original field theory 
[9] . As a final note, we mention that the combination of 
anisotropic dynamics and a conservation law introduces 
power law correlations at all T > T c [3], a manifesta- 
tion of the relaxation of the detailed balance constraint. 
These correlations are revealed by the structure factor, 
which has a discontinuity singularity at the origin. In this 
sense, even the 'disordered' phase is quite non-trivial. 

Turning to multispecies versions, the simplest ('SHZ') 
model [6] allows two different types of particles, distin- 
guished only by their interaction with the external field. 
Positive (negative) particles are biased to hop in the pos- 
itive (negative) y-direction, and interact only via an ex- 
cluded volume constraint. The temperature is absorbed 
into E and the only parameters are E, the overall mass 
density m, and the overall charge density (i.e., the density 
difference of the two species), /. Here, the mechanism for 
ordering is the mutual volume exclusion of the particles, 
so that at sufficiently strong E and large m, the system 
locks into a high density strip perpendicular to E, with 
positive and negative particles blocking each other. For 
non-zero charge density /, this strip is found to drift in 
the direction of the minority species [8] . Depending on 
where the phase boundary is crossed, first-order or con- 
tinuous transitions are observed [14, 15]. Various other 
remarkable properties have been discovered. For a range 
of aspect ratios, configurations with non-zero winding 
number ('barber poles') are quite frequently observed, in 
addition to the usual transverse strips, raising the possi- 
bility of bistability [16]. Power law correlations charac- 
terize the disordered phase, with directionality- dependent 
exponents [17]. Another subtle issue concerns the lower 
critical dimension: While an exact solution for a strictly 
one-dimensional model, characterized by a single 'lane' 
parallel to E, precludes a transition [18], Monte Carlo 
data for a 'two-lane' model indicate the presence of a 
macroscopic cluster in finite systems [20]. Very subtle 
finite-size effects control the decay of this cluster in the 
thermodynamic limit [21-23]. 

In this paper, we consider the two-species model at 
finite T and E, where interparticle interactions are ex- 
pected to play an important role. By varying T, E, 
and /, the fraction of the total population which are of 
the minority species, we can interpolate smoothly from 
the KLS model to the (non-interacting) two-species SHZ 
model. Hence, we expect a competition between the two 
types of ordered configurations - parallel vs transverse 



strips - favored by these two limits. As / varies from 0.0 
(KLS model) to 0.5 (equal numbers of each), there should 
be some critical / where the preferred order switches. 
To explore these phenomena in more detail, we map out 
the phase diagram in E, f and T, for a range of system 
sizes. The energy scale is set by our choice of the inter- 
particle attraction J, and the overall mass density m is 
fixed at 0.5 so that the Ising critical point remains ac- 
cessible. Many questions arise in connection with earlier 
work. How do nearest-neighbor attractions modify the 
two-species transition? What will be the effect of a few 
'impurities' (i.e., minority particles) on the KLS transi- 
tion? At what concentration do the 'impurities' become 
relevant and change the nature of the transition? Prelim- 
inary results, focusing on a restricted parameter space, 
were already reported in [24]; here, we explore a much 
wider parameter range, including several system sizes. 
We will be able, if not to answer these questions fully, 
then to at least suggest the character of their resolution 
sufficiently to guide further research. 

Our main results are as follows. At fixed E and suffi- 
ciently small /, a line of continuous transitions emerges 
from the pure KLS (/ = 0.0) point, in the f-T plane. 
This line separates the disordered phase from an ordered 
one, characterized by a particle-rich strip parallel to E. 
As we increase /, we encounter a bicritical point, where 
the transition line splits into a line of continuous order- 
disorder transitions, from disorder into a strip transverse 
to E, and a line of first-order transitions along which 
transverse and parallel order coexist. If we fix / and 
lower T, we first observe the transition from disorder into 
the transverse strip, followed by a transition into parallel 
order. This topology persists at higher E, except that 
all lines are shifted to slightly higher temperatures. The 
size-dependence of the phase diagram is subtle, since the 
main features are controlled by different scaling variables. 
On the one hand, the transition into the transverse strip 
is controlled by the effective drive L y E/T where L y is the 
system size in the drive direction. On the other hand, the 
bicritical point appears to depend on the scaling variable 
L y f which translates into the number of rows (transverse 
to E) which can be filled with the minority species. Fi- 
nally, the pure KLS point requires finite-size scaling at 
fixed shape factor A oc L y jL\ [10], in two spatial dimen- 
sions. 

The remainder of the paper is organized as follows. We 
first describe in detail the microscopic model and the ob- 
servables which are used to locate the different phases. 
We then present our simulation results, beginning with 
the structure of typical configurations in different parts 
of parameter space and their associated order parame- 
ters. By monitoring the signatures of first and second 
order transitions, we compile a cut through the phase 
diagram at fixed E, with variable / and T. The phase 
boundaries and their dependence on system size are ana- 
lyzed in some detail. To complete the picture, we present 
two cuts at different but fixed temperatures, crossing the 
phase boundaries by varying E and /. We conclude with 
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a brief summary and a discussion of some open questions. 



II. MICROSCOPIC MODEL AND 
OBSERVABLES. 

We consider periodic square lattices of size L x x L y , in 
two spatial dimensions, with E parallel to the (positive) 
y-axis. A configuration is specified by the set of occu- 
pation variables, {a (r)}, where a (r) takes three values, 
±1,0 denoting a positive (negative) particle or a hole at 
lattice site r. Often, we will only need to distinguish par- 
ticles from holes, via n (r) = \a (r)|. All lattices are half 
filled, i.e., m = (i^iy) -1 n ( r ) = 1/2, so that the 
Ising critical point remains accessible. An important pa- 
rameter is the fraction of negative particles (the 'minor- 
ity species') in the system: / = (mL x L y )~ ^2 r ^-i,a(r)- 
Clearly, we only need to consider the sector < / < 0.5, 
from having no negative particles at all to equal num- 
bers of each species. For later reference, we also intro- 
duce the charge density, q = (L x L y )~ x ^ r a (r) = m — f. 
The nearest-neighbor attraction is modelled by the Ising 
Hamiltonian 



H = -AJ n(r)n(r') 



(1) 



<r,r') 



We choose attractive interactions, J > 0, regardless 
of species. While many other choices are possible and 
interesting, ours provides maximum linkage to known 
cases: Ising, KLS and SHZ. The Monte Carlo dynam- 
ics conserves the number of each species and is speci- 
fied as follows. An update attempt begins by picking a 
bond at random. If the bond connects a particlc-holc 
pair, the contents are exchanged with the Metropolis 
rate min{l,exp [- (AH - SyEa (r)) /T]} [25]. Hence, at 
E = we recover the equilibrium Ising model with con- 
served magnetization, coupled to a heat bath at temper- 
ature T. We will set J = 1 and measure E in units of 
this (arbitrary) energy scale. T will be quoted in units of 
the Onsager temperature, T C (E = 0). The change in in- 
coordinate, due to the proposed move, is denoted by Sy, 
and AH is the associated change in internal energy. The 
term SyEa (r) models the gain or loss of energy from the 
coupling to E; if Sya (r) is positive (negative) the move 
is favored (unfavored). Our model, in which E and T 
are varied independently, raises an interesting issue. If 
the ratio E/T is quite large, it becomes almost impos- 
sible for particles to hop backwards. In a finite system, 
this implies that a relatively small fraction of the mi- 
nority species - provided the 'right' fluctuation occurs 
- is sufficient to form a stable blockage. Even though 
such a fluctuation becomes less probable in a larger sys- 
tem, the dynamics nevertheless becomes nonergodic in 
the limit E/T — > oo. In principle, this can be avoided 
by introducing, e.g., a small probability for particles to 
exchange places [15]. To limit the number of parameters, 
we circumvent these problems here by considering differ- 
ent initial configurations and a range of system sizes. 



The dynamics is diffusive, and therefore conserves both 
charge and mass density. Though the local effect of the 
external field is analogous to the effect of an electrostatic 
potential on electric charges, the boundary conditions ex- 
clude the possibility of a global Hamiltonian description. 

As overall density is conserved, we expect ordered 
configurations to be strips of higher density coexisting 
with strips of lower density. We therefore introduce the 
Fourier transform of the local mass variable, 



s(m x ,m y ) 



L x L y — ' 



n(x,y)e 



2iri(m x x/L x +m y y/Ly) 



(2) 



which is labelled by (integer) wavenumbers, m x — 
0, 1, L x , m y = 0, 1, L y — 1. The structure factor, 



S(m x ,m y ) = (\s(m x ,m y 



(3) 



then serves as a good order parameter, since it is sensi- 
tive to mass-segregated strip configurations. For example 
5(1,0) will be O (1) for a strip aligned with the field, 
characteristic of KLS order; similarly, S (0, 1) will detect 
a strip transverse to E which develops in the SHZ (two- 
species) model; and both are normalized to O (l/L x L y ) 
for a disordered configuration. We also monitor a 'sus- 
ceptibility', i.e. the fluctuations of the order parameter: 

A(m x ,m y ) = L x L y (^\s(m Xl m y )\ 4 ^ (4) 

.. . 21 

\s(m x ,m y )\ J 

We note that S(m x ,m y ) involves the Fourier transform 
of the mass variable and is therefore not sensitive to 
any charge-segregated structures. Replacing n(x, y) by 
cr(x,y) in Eq. (3) generates structure factors which re- 
spond to charge inhomogeneities. We have monitored 
these and their fluctuations throughout, and found that 
their behavior is consistent with the mass-based quanti- 
ties. 

When S is calculated, the average is taken over mul- 
tiple steady-state configurations of a Monte Carlo run, 
with a typical run lasting 0.8M (8 x 10 5 ) Monte Carlo 
steps (MCS), and 2L x L y bond update attempts per 
MCS. Data are collected every 400 MCS; fluctuations 
of observables indicate that this interval is sufficient to 
produce uncorrelated data in the largest (60 x 80) sys- 
tems considered. Typically, the initial 0.2M MCS are 
discarded to ensure that data are taken from the steady 
state. Near critical points and at low temperatures these 
numbers require modification, due to long correlation 
times and long-lived metastable states. In such cases the 
only recourse is a careful analysis of individual, very long 
runs. When that is necessary we will measure a quantity 
closely related to S: 



s(m x ,m y ) = \s(m x ,m y )\* 



(5) 



which measures the type of order present in a single con- 
figuration. We can then track s for different m x s and 
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FIG. 1: Phase diagram in / and T for E = 2.0. Triangles and 
diamonds are boundaries in the 40 x 40 system, x's are for 
the 60 x 60 system. 

m y 's over the course of a run and see precisely how the 
averages are generated. 

Finally, at each bond update we tally the quantity 
Sya (x, y), which is then averaged over a run to give the 
charge current j. However, it is not particularly illumi- 
nating to compare j for different values of the effective 
drive E/T. We are more interested in a quantity which 
is a property of the gross structure of the steady state, 
namely, the conductivity k defined via j — kE/T. 

Now that we have described the various quantities 
which will be used to probe the behavior of our model, 
we turn to the presentation of the data. 



III. RESULTS 
A. Phase diagram in / and T 

In this section we seek the location and character 
of transitions by scanning in / and T at fixed drive 
E. We choose E = 2.0 since this intermediate value 
still allows for a significant fraction of backward jumps, 
thus avoiding the spurious metastable configurations dis- 
cussed above. At the same time, it is large enough to 
induce measurable currents and other clear signatures of 
far-from-equilibrium behavior. The two order parame- 
ters S (1, 0) and S (0, 1) and their fluctuations are mon- 
itored in order to identify the different phases. Large 
peaks in their fluctuations, or the presence of hysteresis, 
are used as indicators of continuous vs first order tran- 
sitions, respectively. For clarity, we first present a quick 
overview of the topology of the phase diagram, and then 
turn to the details of the data which underlie this picture. 

Fig. 1 shows the phase diagram in the f-T plane, at 
E = 2.0, for two different system sizes. Three phases are 
found: a homogeneous, disordered phase (DO), a trans- 
verse strip (TS) phase as in the two-species model, and 
a parallel strip (PS) phase as in the KLS model. The 
value of / determines which phase is observed: at / = 




FIG. 2: Configurations of the 40 x 40 lattice at / = 0.50 for 
four temperatures. Upper left, T = 6.0; upper right, T = 
3.25; lower left, T = 2.64; lower right, T = 1.0. Positive 
(negative) particles are white (black); E points upwards. 

there is only one species of particles, and a single transi- 
tion is observed from disorder into the parallel strip. As 
/ increases, this transition persists until the number of 
the minority species is sufficient to create a blockage and 
form the transverse strip. From here on, two transitions 
are observed: from disorder into the transverse strip, and 
at a lower temperature from the transverse strip into the 
parallel strip. Upon increasing / further, only the DO- 
TS transition can be detected. Although the TS-PS line 
cannot have a critical endpoint for reasons of symmetry, 
for / > 0.10 it occurs at such a low temperature that it 
cannot be observed in simulations of a reasonable length. 

Now that we have briefly discussed the phase diagram 
we can look in more detail at the phases and their bound- 
aries. Before presenting the data (structure factors, their 
fluctuations, currents) we will begin with some pictures 
of typical configurations at various points in the phase di- 
agram. This way the reader can develop some intuition 
about the model. All pictures show 40 x 40 lattices, with 
E pointing up. White (black) pixels are positively (neg- 
atively) biased particles, blue-green (gray in print) pixels 
are holes. We will begin at the right side of the phase di- 
agram and explore how typical configurations change as 
we move left, decreasing the fraction of the population 
which belongs to the minority species. 

In Fig. 2 we present four configurations at / = 0.50: 
The first picture is at T = 6.0, well above T c = 3.3. 
Unsurprisingly, this configuration lacks any visible struc- 
ture since it falls deep into the homogeneous phase. This 
will be the only picture presented for the homogeneous 
phase, as the only visible feature which changes with / 
is the ratio of white to black pixels. The next picture 
is at T — 3.25, and now we begin to see the two-species 
type phase separation. Lowering the temperature further 
to T — 2.64, the horizontal strip appears quite clearly, 
though it remains diffuse at the boundaries, and there 
are many 'travelers' (isolated particles) moving through 
the remainder of the system. At this temperature back- 
ward hops of particles are not too improbable, occuring 
with a rate exp(~E/T) = 0.47, and thus a fair number of 
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FIG. 3: Configurations of the 40 x 40 lattice at / = 0.075 
for four temperatures. Upper left, T = 1.77; upper right, 
T = 0.95; lower left, T = 0.84; lower right, T = 0.78. 

holes are able to enter the strip and allow particles to slip 
through the blockage. Finally at T = 1.0 we are deeply 
inside the ordered phase, and the boundaries of the strip 
are very sharp and travelers are few. At this tempera- 
ture, holes are rarely able to penetrate into the interior 
of the strip. The irregular shape of the interface between 
the two species is a result of the quench from disorder; 
a more gradual lowering of the temperature would result 
in a smoother interface. 

Lowering / to 0.075 (Fig. 3), only 1.5 rows of the mi- 
nority species are left. The first frame shows a config- 
uration just below criticality at T = 1.77. Though a 
blockage can still form, the strip is no longer symmet- 
ric with respect to + and — . This leads to a drifting 
of the strip: Occasionally, the rather thin blockage of 
the minority species is opened by backward hops, and 
the majority species pours through. These particles then 
travel quite rapidly around the periodic lattice and at- 
tach to the back of the majority blockage, the net result 
being an upward drift of the strip. Lowering T further to 
0.84 (third frame) it appears that interfaces parallel to 
E are becoming favorable; this type of configuration is 
common at these intermediate values of /: here, parallel 
and transverse strips compete with each other. Indeed 
the final frame (T = 0.78) shows the preferred low tem- 
perature configuration: a single strip of mixed charge 
parallel to E, raising the possibility of a sequence of two 
transitions as a function of T. 

Now that we have developed a qualitative notion of the 
various regions in the phase diagram, we turn to a more 
quantitative analysis of the phases and their boundaries. 
As in the preceding section, we present the order param- 
eters and their fluctuations as a function of T, for a range 
of /, though with one important difference: we will show 
data for at least two system sizes. While this cannot re- 
place a high-precision finite-size scaling analysis of the 
transitions, it is intended to provide a rough picture of 
how the fluctuations of the order parameter scale with 
system size. 

We begin as before at / = 0.50, with equal numbers 
of each type of particle. In Figs. 4 and 5, 5(0, 1) and its 




FIG. 4: S (0, 1) as a function of T for / = 0.50. Open (filled) 
triangles are for the 40 x 40 (60 x 60) system. 
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FIG. 5: A (0, 1) as a function of T for / = 0.50. Open (filled) 
triangles are for the 40 x 40 (60 x 60) system. 

fluctuations are plotted as a function of T for two system 
sizes (40 x 40 and 60 x 60); 5(1, 0) is not shown as the only 
transition here is from the homogeneous phase into the 
transverse strip. In both systems, 5(0, 1) goes smoothly 
to zero as T is increased. A clean peak in A(0, 1) is also 
observed in each system, increasing in amplitude with 
the system size. These two observations are consistent 
with a continuous transition into the transverse strip at 
/ = 0.50, and we therefore use the location of the peak 
in A(0, 1) to locate the phase boundary in Fig. 1, with 
T c (L = 40) = 3.35 and T c (L = 60) = 4.64. 

Reducing / to just a few rows of the minority species 
introduces some new features into the data, which we 
describe briefly before continuing our analysis. Here, 
we still observe the transition described above, with the 
same behavior of 5(0, 1) and A(0, 1). However, now that 
the strip is no longer symmetric, it drifts. The ordered 
phase therefore fluctuates significantly, and A(0, 1) de- 
velops a shoulder at about 1 /3 of the value it reaches at 
the transition. We will see below that the fluctuations of 
the ordered phase can obscure the transition in smaller 
systems. 
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Upon reducing / further we observe a sequence of two 
transitions as a function of T . Fig. 6 shows both order 
parameters, 5(0,1) and 5(1,0) for 1.5 rows of the mi- 
nority species: / = 0.075 in the 40 x 40 system. Also 
shown in Fig. 7 is A(0, 1) for 1.5 rows of the minority 
species in both the 40 x 40 and the 60 x 60 system. Wc 
have omitted the 60 x 60 data for the order parameters to 
keep the plot uncluttered. As before, 5(0, 1) and A(0, 1) 
signal a continuous transition into the TS phase as T is 
lowered, though the signal in A(0, 1) is much more pro- 
nounced in the 60 x 60 system. And also as before, there 
are significant fluctuations associated with this phase, 
due to strip drifting. Specifically, in the 60 x 60 sys- 
tem A(0, 1) actually has a broad secondary peak in the 
ordered phase. The large magnitude of this signal is is 
quite unexpected and awaits a satisfactory explanation. 
For now, we only note that lowering T increases the ef- 
fective bias, E/T , and therefore enhances fluctuations 
associated with the drive. Lowering T further we ob- 
serve 5(0, 1) falling abruptly, while 5(1, 0) climbs rapidly, 
suggesting a discontinuous transition from the horizontal 
strip (TS) into the vertical strip (PS). In the neighbor- 
hood of such a transition, one expects to see metastabil- 
ity of the unfavored phase, and this is indeed the case as 
shown in Fig. 8. Here we have plotted time traces (as op- 
posed to configurational averages) of structure factors for 
individual configurations, s(l,0) and s (0, 1), defined in 
Eqn. (5). Clearly, when s (1, 0) (s (0, 1)) = 1 the configu- 
ration is a perfect vertical (horizontal) strip. Sufficiently 
close to the transition, the time traces reveal the expected 
behavior, as the system switches between the two ordered 
phases. Notice the length of the run shown: 40M MCS, 
which is a factor of 40 longer than typical runs, indi- 
cating that the lifetimes of metastable configurations are 
already quite long even in the 40 x 40 system, rendering 
such behavior nearly unobservable in the 60 x 60 system. 

At smaller values of / we are nearing the junction of 
the three phase boundaries, which considerably compli- 
cates the analysis of data from small systems in a couple 
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FIG. 6: S (0, 1) and S (1, 0) as a function of T for / = 0.075 
in the 40 x 40 system. Triangles (squares) are for S (0, 1) 
(S (1, 0)). Lines and dashes are provided to guide the eye. 
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FIG. 7: A (0, 1) as a function of T for / = 0.075. Open (filled) 
triangles are for the 40 x 40 (60 x 60)system. 



of ways. The sequence of transitions (DO-TS followed 
by TS-PS) becomes difficult to resolve, as they are quite 
close in temperature, and massive fluctuations from the 
first-order TS-PS line may wash out the signal in A(0, 1) 
which locates the continuous DO-TS line. And if the 
junction of the three lines is indeed a nonequilibrium bi- 
critical point, we can expect unexplored finite-size effects 
to interfere with the analysis. We can, however, make 
some progress based on the assumption that the relevant 
control parameter near the bicritical point is the number 
of rows of the minority species. This hypothesis will be 
treated in more detail below, in the sections on scaling 
arguments. 

At precisely one row of the minority species it is no 
longer possible to accurately resolve the two transitions 
in the 40 x 40 system. A weak transverse ordering is 
observed, with 5(0,1) reaching at most 40% of perfect 
order. In the vicinity of the PS-TS first-order transi- 
tion, huge fluctuations associated with switching between 
metastable configurations are observed, which wash out 
the signal of the DO-PS transition. However, it is inter- 
esting that these transitions can be resolved in larger sys- 
tems at precisely one row of the minority species: There, 
these two transitions are sufficiently far apart in temper- 
ature since (as we will see below) T c increases with L y 
across the DO-TS transition. This likely explains why 
the DO-TS transition is not observed in the 40 x 40 sys- 
tem. 

Below one row of the minority species we no longer ob- 
serve the transverse order, though we caution that this 
may be strictly correct only for the finite system. With / 
just below a single row of the minority species, the high 
temperature phase is homogeneous and the low temper- 
ature phase is the parallel strip. In the vicinity of the 
transition, huge fluctuations are observed in 5 (0, 1) and 
especially 5 (1, 0), where the fluctuations are an order of 
magnitude larger than the signal at the DO-TS bound- 
ary. In this region neither 5(1,0) nor 5(0,1) posesses 
a well-defined average; timetraces indicate vigorous com- 
petition between the two ordered phases. We conjecture 
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FIG. 8: Timetrace at T = 0.832. Time in units of 0.2M MCS is plotted on the horizontal axis. The values of s (0, 1) and s (1, 0) 
(triangles and squares, respectively) are plotted on the vertical axis. 



that we are close to the bicritical point in the finite sys- 
tem, and are therefore unable to resolve the transition 
without some knowledge of scaling to guide the analy- 
sis. At smaller / we are farther from the bicritical point, 
and the complications from the presence of the minority 
species are less severe. Fig. 9 shows 5 (1, 0) for two sys- 
tem sizes at exactly one-half row of the minority species. 
We have not shown 5 (0, 1), as the signal has become in- 
significant. 5(1,0) shows the low-temperature configu- 
ration to be a single parallel strip, the smooth approach 
to zero again suggesting a continuous transition. And 
indeed, the data for A (1, 0) is consistent with this con- 
jecture, with a sharp peak which increases with system 
size. The amplitude of the peak is smaller by a factor 
of five than the peaks observed near the bicritical point, 
suggesting that it is associated with fluctuations about a 
well-defined average, as opposed to transitions into and 
out of the ordered phase. Timetraces support this con- 
jecture. 

Charge currents were also measured in order to look 
for signatures of the various phases and transitions. We 
note however that as we are changing T, we are also 
changing the effective bias e = E/T, and therefore a 
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FIG. 9: 5(1,0) as a function of T for / = 0.025. 40 x 40 
(60 x 60) data are shown by open (filled) squares. 



more appropriate quantity is the conductivity, k = j/e. 
In Fig. 10 we plot function of e for each / dis- 

cussed in this and the following sections; we present only 
data for the 40 x 40 system as there are no significant 
differences between the two system sizes. At tempera- 
tures below T = 1.4 (e > 1.4) the conductivity vanishes 
in the / = 0.50 systems; as the temperature is raised 
the effective drive is reduced, and backward hops occa- 
sionally occur, allowing a small current to trickle through 
the blockage. Upon raising the temperature further, the 
conductivity changes slope at T = 3.6 in the 40 x 40 
system, which does not coincide with the phase tran- 
sition. Rather the maximal conductivity occurs in the 
disordered phase, but at a temperature which is not too 
large, so that backward hops are not so common as to 
begin reducing the current. The transition apparently 
corresponds instead to the inflection point (T = 3.4), 
where the curvature changes sign. At / = 0.10, the con- 
ductivity has a slope discontinuity in the 40 x 40 sys- 
tem at T = 2.6, slightly above the critical temperature 
T c = 2.4. Though there should be an inflection point in 
k near the transition, our data are not precise enough to 
locate it. At f — 0.075 and below, the conductivity drops 
smoothly to zero with T, showing no indication of either 
transition. Apparently, 5 and A are far more sensitive 
to the transitions in our model. 

We close this section on the f-T phase diagram with 
a summary of the results. The picture at higher / is 
clear: a clean continuous transition into the horizontal 
strip, with T c decreasing with /. When / is reduced 
to approximately three rows of the minority species, the 
signal of the transition remains clear, though it now sits 
atop a shoulder of fluctuations of the ordered phase. At 
yet smaller /, a second transition appears between the 
two ordered phases at lower T; it has the characteristics 
of a first-order transition. At even lower /, at approxi- 
mately one row of the minority species, both fluctuations 
perpendicular and parallel to E become so violent that 
the DO-TS transition is only seen in larger systems, as 
the two transitions nearly overlap in the 40 x 40 system. 
Close to the / = point, the transition is once again 
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clean and apparently continuous into a vertical strip of 
mixed charge. Of course, all these statements are based 
on an analysis of finite systems. In order to draw robust 
conclusions, a more systematic analysis of larger samples 
is required. 



B. Phase diagram in / and E 

In the preceding section we studied a slice of the phase 
diagram at constant E, varying the fraction of the minor- 
ity species and the temperature. Varying T effectively 
varies both the strength of particle-particle attractions 
and the strength of the bias, since the relevant quanti- 
ties in the rates are J/T and E/T. In this section we 
consider a different cut through the phase diagram. By 
varying E and / at fixed T, the interparticle attractions 
are held constant while the strength of the bias is var- 
ied. In this way we can study directly the competition 
between the drive and the attractive interactions. In the 
following, we choose a value for T, and then scan in E 
for several values of /. The temperatures are chosen by 
reference to the KLS temperature: T = 2.0 is above the 
critical temperature of the KLS model at saturation, and 
T = 1.2 is at the critical temperature for E = 2.0, stud- 
ied in the previous section. At this stage, we have only 
data for 40 x 40 systems, and are therefore as yet un- 
able to speculate on results for larger systems. However, 
they cast a new light on the more detailed results of the 
previous sections. 



1. T = 2.0 

As before, we first survey the phase diagram with the 
help of some typical configurations. The qualitative pic- 
ture will then be made more quantitative in the next 
section, by examining the behavior of currents and order 
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FIG. 10: Conductivity as a function of effective drive E/T for 
several /. / = 0.50, diamonds; / = 0.10, squares; / = 0.075, 
triangles; / = 0.05, circles; / = 0.04, +'s. Lines have been 
added to the largest and smallest systems to guide the eye. 




FIG. 11: Configurations for several / at E = 20.0, T = 2.0. 
Upper left, / = 0.50; upper middle, / = 0.10; upper right, 
/ = 0.075; lower left, / = 0.05; lower middle, / = 0.04; lower 
right, / = 0.025. 

parameters. 

Fig. 11 shows a series of configurations at various / for 
E = 20.0. The first frame clearly shows the transverse 
strip at / = 0.50, and the absence of travelers suggests 
that the strip is stationary. In the next frame we have 
reduced / to 0.10, reducing the thickness of the minority 
species to exactly two rows. Consequently we now see 
some travelers trickling through a break in the blockage. 
Watching an animation in this region of the phase dia- 
gram reveals an interesting behavior: the strip is mostly 
quiescent, except for a few particles hopping back and 
forth at the particle-hole interface. Aside from the dif- 
ferent ratio of + to — , these configurations look similar 
to the / = 0.50 strip. Then, a sudden large fluctuation 
opens up a hole in the minority blockage: the + par- 
ticles pour through, and the strip fluctuates and drifts 
partway around the lattice, until the blockage is reestab- 
lished. Reducing / further to 0.075 (third frame) we see 
a strip in the middle of one of these fluctuation events. 
In contrast to / = 0.10 where such large fluctuations are 
relatively rare, the situation is now reversed; i.e., the qui- 
escent periods become less frequent. In the fourth frame, 
we set / — 0.05, and while the strip is still clearly visible 
it now drifts continuously. The final two frames show 
/ = 0.04 and 0.025. Now there is no longer any clear 
evidence of phase separation. This rough picture is con- 
sistent with our earlier investigation of the two-species 
transition, where we observed transverse order at and 
above a single row of the minority species. 

Fig. 12 presents the phase diagram at T = 2.0. The 
boundary separates a horizontal strip at high E and / 
from a homogeneous phase at small E and /. As we are 
above the critical temperature for the KLS model at sat- 
uration bias, the vertical strip does not appear at any E; 
at low / (where we might otherwise expect to see such 
ordered configurations) the system simply remains dis- 
ordered for any E and /. Figs. 13 and 14 show S (0, 1) 
and A (0, 1) for several values of /; in the interest of 
clarity A (0, 1) is plotted for only four /'s. As E is in- 
creased the system orders into a transverse strip, with 
S (0, 1) saturating at smaller values as / is decreased. At 
/ = 0.05 (exactly one row of the minority species) S sat- 
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FIG. 12: Phase diagram in / and E/T at T = 2.0. Disorder 
(DO) is observed at small / and E, the transverse strip (TS) 
dominates at high / and E. The error bars are smaller than 
the size of the data points, unless indicated. 
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FIG. 13: S (0, 1) as a function of E/T for several /. / = 0.50, 
diamonds; / = 0.10, squares; / = 0.075, triangles; / = 0.05, 
circles; / = 0.04, +'s; / = 0.025, x's. 
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FIG. 14: A (0, 1) as a function of E/T for several /. The 
symbols are the same as in the previous figure. A line has 
been added to the / = 0.075 data to guide the eye. 



urates at only 0.28, indicating that the transverse order- 
ing is rather weak, though comparison with the data for 
/ = 0.025 shows dramatically different behavior. Here 
S (0, 1) reaches a maximum of only 0.01, and the behav- 
ior can hardly be called 'saturation'. A (0, 1) also signals 
a transition, though the clean, sharp spike at / = 0.50 
becomes a broad bump at / = 0.05, and shows no sig- 
nal at / = 0.025. The susceptibility also indicates a 
difference in the ordered phases at different /: at large 
/ increasing E suppresses fluctuations, while at smaller 
/ (when the strip begins drifting) increasing E enhances 
fluctuations. It is important to note that the fluctuations 
at high E are fluctuations about the ordered phase, as 
A (0, 1) is always 2 to 3 orders of magnitude larger than 
A (1,0). Though the data is not included in the plots, 
we have checked the behavior of the ordered phase for E 
as high as 40. The fluctuations for small / (/ < 0.10) 
saturate and bounce around a well-defined average, while 
for larger / they are suppressed. This is true for / = 0.05 
and greater; at / = 0.025 the magnitude of fluctuations 
in either direction are comparable. We also examined 
the conductivity across the phase boundary. We do not 
present the data here, as it is qualitatively similar to the 
data in Fig. 10. Again, the transition apparently coin- 
cides with the inflection point of the conductivity. 



2. T = 1.2 

In the preceding section we found a simple phase dia- 
gram, with a single boundary separating two-species or- 
der from the homogeneous phase. Since the tempera- 
ture was chosen well above the maximum KLS transition 
temperature, Tkls {E = oo) = 1.414, the KLS ordered 
phase could not be observed. Now, we lower the temper- 
ature to T = 1.2 < Tkls an d explore the corresponding 
(f-E) slice of the phase diagram (Fig. 15). As long as 
we remain at / > 0.10, we observe a transition similar 
to the one at T — 2.0: from the homogeneous phase into 
two-species order. In contrast, for / < 0.04 the KLS 
transition is observed, since the minority species is too 
scarce to form a blockage and T < Tkls (E = oo). Be- 
tween these two limiting values of / we are again in the 
vicinity of the bicritical point, and the situation becomes 
complicated due to the competing types of order. 

Fig. 16 shows typical configurations at two different 
values of / for various E. For / > 0.075 these configu- 
rations look much as they did at T = 2.0, so they need 
not be included. The / = 0.075, E = 20.0 configura- 
tion shows some very interesting structure, almost 'equal 
parts' KLS and two-species order, suggesting some seri- 
ous competition between the two phases. We stress that 
this is a typical configuration. When / is reduced to a 
single row (/ = 0.05) this competition is reduced, and we 
see instead a KLS phase with some local two-species or- 
der. This trend continues upon reducing / to zero, where 
at high E the KLS order is observed. The other panels 
show / = 0.075 and 0.05 at smaller values of E, E = 4.0 
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and 2.4 respectively. These E values were chosen be- 
cause they maximize the two-species order for these f's. 
In each case the strip drifts rapidly around the lattice. 
Interestingly, as the majority species is piled onto the 
back of the drifting strip it builds long fingers, leading 
to a very irregular interface. Pictures at smaller /'s are 
not shown, as they do not form the transverse strip at 
smaller E. In the next section we locate these boundaries 
and study the phases in greater detail by examining the 
behavior of S, A, and k. 

In Fig. 17 we plot S for the high-/ phases. For / = 0.50 
and 0.10, S (0, 1) shows the system ordering into the two- 
species phase much as in the previous section. There is, 
however, one notable new feature for / = 0.10: we see a 
slight suppression of S (0, 1) over an intermediate range 
of E, from about E — 6 to E — 20. In this range back- 
ward hops occur frequently enough to occasionally open 
a hole in the minority blockage, allowing the majority 
species to pour through until the blockage reforms. In- 
creasing E reduces the likelihood of such events, until 
the strip is perfect except for some surface diffusion at 
the particle-hole interface. At / = 0.075 the behavior 
changes dramatically. After maximizing the two-species 
order at E = 3.0, increasing E further suppresses S (0, 1) 
and enhances S (1, 0), until both saturate below 0.2. The 
system can hardly be said to have switched to KLS or- 
der, though neither is it precisely disordered either. To 
understand this peculiar behavior, we recall that T (and 
therefore J/T) is held constant, so that increasing E fa- 
vors interfaces parallel to E. Yet, the number of mi- 
nority particles is large enough to form bubbles of local 
two-species order, as one can discern from the top left 
panel of Fig. 16. Apparently the competition between 
the two phases is very balanced in this small system. It 
would be interesting to simulate larger systems and ex- 
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FIG. 15: Phase diagram in E/T and / for T — 1.2. A con- 
tinuous boundary separates DO from TS (diamonds) as well 
as DO from PS (triangles) . x 's indicate a possible boundary 
between TS and PS. A few typical error bars are shown. 
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FIG. 16: Typical configurations for the T = 1.2 plane. Upper 
left, / = 0.075, E = 20.0; upper right, / = 0.05, E = 20.0; 
lower left, / = 0.075, E = 4.0; lower right, / = 0.05, E = 2.4. 



plore whether this type of 'phase competition' persists, 
or whether the KLS order eventually becomes stable. 

Additional information is provided by the susceptibil- 
ities. Observations of A (0, 1) are consistent with a con- 
tinuous transition from disorder into the transverse strip 
at small E. Seeking a signature of the KLS phase, we 
show A (1, 0) and A (0, 1) for / = 0.075 in Fig. 18. Here 
we observe a second peak in A (0, 1) at E = 7.2, which 
corresponds to S (0, 1) ~ 0.4 in Fig. 17. This peak is 
rather broad and its amplitude is more than twice that 
of the first peak (associated with the DO-TS transition). 
Perhaps it suggests a first-order transition which would 
be observed in a larger system, separating the two-species 
phase from the KLS phase. Finally, at higher E both 
A (0, 1) and A (1, 0) fluctuate around nonzero values, re- 
flecting the fluctuations of the competing phases. We 
note no signal of a transition in A (1, 0). 

We now continue by looking at the data for / = 0.05, 
shown in Fig. 19, where we have plotted S for each type 
of order. We observe that the transverse strip only makes 
a brief appearance, with S (0, 1) reaching a maximum of 
only 0.35. At higher E, the parallel strip stabilizes and 
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FIG. 17: 5" (0, 1) as a function of E/T for several /. / = 0.50, 
diamonds; / = 0.10, squares; / = 0.075, filled triangles. Open 
triangles indicate S (1, 0) for / = 0.075. 
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reaches a value slightly above 0.4. However, we observe 
no clear signature of a distinct transition between the 
two types of order. There is a vague remnant in A (0, 1) 
of the double peak structure seen at / = 0.075, and 
then a broad, messy bump in A (1,0) before it settles 
down to capture the fluctuations of the ordered phase. 
Only further investigation of larger systems can deter- 
mine whether these weak signals in fact indicate a tran- 
sition. It is clear, however, that the high E phase here is 
not the 'phase competition' seen previously at / = 0.075, 
since neither 5 (0, 1) nor A (0, 1) contain any trace of it. 
Upon reducing / below 0.05 we find that the transverse 
strip has essentially disappeared. Meanwhile, 5(1,0) 
looks similar to Fig. 19, saturating to about 0.50, indi- 
cating that the KLS strip has formed. Apparently the 
system is unable to completely order at any E, since 
Tkls (E = oo) = 1.4. We stress that this behavior is due 
to the proximity of the KLS phase transition, not due to 
some residual competition with the transverse strip. Fi- 
nally, we note that the behavior of the susceptibility is 
consistent with a continuous transition into the parallel 
strip. 

The conductivities are shown in Fig. 20. At / = 0.50 
and 0.10, k has much the same form as seen before, de- 
creasing monotonically with E. This is no longer the 
case when / is reduced to 0.075. Now k develops a broad 
minimum which coincides with the appearance of the 
two-species order, followed by a small peak at E = 8 
before falling off with E. The peak coincides with the 
crossing of 5(0,1) and 5(1,0), which suggests that the 
current is maximized by the 'phase competition'. With 
the current maximized, increasing E simply reduces k; 
indeed a glance at the raw data for the current shows 
that this is the case. The data for / = 0.05 shows the 
same nonmonotonic structure as for the / = 0.075 data. 
Again, the peak corresponds to the crossing of the two 
order parameters, and then falls off with increasing E. 
This similarity suggests that perhaps the high E phase 
at / = 0.075 is in fact the parallel strip, and such behav- 
ior would be observed in a larger system. For values of 
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FIG. 19: S(0, 1) (filled circles) and 5(1,0) (open circles) as 
a function of E/T for / = 0.05. 



/ below 0.05 we no longer observe the secondary peak, 
nor do we observe any phase competition. We must, 
however, leave questions about the nature of the high-i? 
phase open to further study. 



C. Scaling Arguments 

Now that we have surveyed the phase diagram in some 
detail, we turn to look closely at some of the boundaries 
with the help of some scaling arguments. Considering 
Fig. 1 again, we note two potentially troubling features. 
First, there is a shift in T c across the DO-TS boundary of 
about 50% between the two system sizes. We will charac- 
terize this shift using a mean-field scaling argument. Sec- 
ond, we note that the bicritical point has shifted towards 
the / = axis in the larger system. It has been alluded 
to before that the number of rows of the minority species, 
rather than the fraction /, might be the controlling vari- 
able. We will investigate this suggestion more carefully 
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FIG. 18: A (0, 1) (filled triangles) and A (1, 0) (open triangles) 
as a function of E/T for / = 0.075. 
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FIG. 20: Conductivity in the T = 1.2 plane as a function of 
E/T for several f.f = 0.50, diamonds; / = 0.10, squares; 
/ = 0.075, triangles; / = 0.05, circles; / = 0.04, +'s. 
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FIG. 21: 5(0, 1) as a function of effective drive for J = 1. 
System sizes: 40 x 40, open triangles; 40 x 60, circles; 40 x 80, 
x 's; 60 x 60, filled triangles; 60 x 80, squares. Lines have been 
added to the largest and smallest systems to guide the eye. 



by considering some larger systems and rectangular ge- 
ometries. Finally we study the / = phase transition, 
using scaling arguments developed for the KLS model. 

The shift in T c with system size is most pronounced at 
/ = 0.50. Previous work on the two-species model with 
J = treated the ordered phase in a mean-field approxi- 
mation by solving equations of motion for the two differ- 
ent charge densities [8, 14]. It was found that the scaling 
functions depend on the combination EL y /T, indicat- 
ing that the effective bias E/T introduces a new length 
scale. This scaling implies an infinite-volume limit in 
which E/T — > as L y — > oo, while keeping EL y /T fixed. 
Earlier analyses of the ordered phase based on these ideas 
have worked quite well, so that we now attempt to ex- 
tend this approach to analyze quantities near criticality 
and for J ^ 0. There is no reason to expect success a 
priori, as both critical fluctuations and nonzero J may 
modify the scaling variables and the mean-field expo- 
nents. In Figs. 21 and 22 we have plotted 5(0,1) and 
A (0, 1) for J = 1.0. (We have divided A (0, 1), Eqn. (4), 
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FIG. 22: A (0, 1) as a function of effective drive for J = 1. 
The symbols are the same as in the previous figure. 



by the volume in order to compare different system sizes 
more easily.) Rather than crossing the phase boundary 
by varying T we have opted instead to vary E since this 
allows us to vary the effective bias E/T at constant in- 
teraction strength, J/T . While the collapse of S (0, 1) 
in Fig. 21 is not perfect, the mean-field scaling argument 
accounts for most of the shift in T c . There is a shift in the 
peak of A (0, 1) of about 3.6% between the largest and 
smallest systems; if the same data were plotted without 
rescaling, the shift in T c is about 52%. Also of note is the 
extremely weak dependence on the transverse dimension, 
as predicted by the scaling argument. We also examined 
the same quantities for J = 0. Interestingly, the data col- 
lapse for J = 1.0 is better than in the J = case. This 
is somewhat surprising, as the original scaling argument 
was derived for J = 0, and we would expect the interac- 
tions to perhaps modify it. Whatever the resolution of 
this puzzle, it is apparent that the scaling argument pre- 
sented here explains the pronounced shift in T c seen in 
the phase diagram: increasing L y requires a correspond- 
ing decrease in the effective bias E/T. 

Another issue concerns the location of the junction of 
the three phase boundaries, shown in Fig. 1. It is clear 
that the junction moves toward the / = axis as the 
system size is increased. In fact, in both systems the 
boundaries merge just below a single row of the minor- 
ity species, which naturally corresponds to a smaller / 
in the larger system. In Fig. 23 we have replotted the 
data from Fig. 1, replacing / with fL y /2, which is sim- 
ply the number of rows of the minority species. Near the 
junction of the three lines we have also included results 
from a few other system sizes with rectangular geome- 
tries. Plotted vs. fLy/2, the junctions of the boundaries 
coincide, within the error bars, for all system sizes, sug- 
gesting that the onset of the two-species order occurs, at 
least in relatively small finite systems, when there are suf- 
ficient minority particles to form a single row. The cru- 
cial question concerns the extrapolation of this result to 
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FIG. 23: The phase diagram in the f-T plane, with / rescaled 
to represent the number of rows of the minority species. Sys- 
tem sizes: filled triangles, 40 x 40; x's, 60 x 60; squares, 
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an appropriate thermodynamic limit. If the system size 
goes to infinity in the most naive way, i.e., L x , L y — > oo at 
fixed aspect ratio L x /L y , the particle density associated 
with a 'single row' vanishes. It is possible that the DO-PS 
transition exists in an infinite volume only at / = 0, and 
any finite density of 'disorder' (i.e., the minority species) 
induces the two-species order. Preliminary studies [26] 
indicate that the minority species does indeed constitute 
a relevant perturbation to the KLS fixed point. We will 
have to leave discussion of this issue to future work and 
for now limit ourselves to statements about finite sys- 
tems. 

At / = there is only one species, and we observe the 
KLS transition at finite E. Though a great deal of study 
has been devoted to this transition at infinite E, there 
has been no detailed work at finite E. Here we present a 
basic finite-size scaling (FSS) analysis of this transition, 
in order to locate T c (E = 2.0), and also to demonstrate 
the subtleties which can arise when studying phase tran- 
sitions with anisotropic, noncquilibrium dynamics. 

Field-theoretic studies of the KLS model [12] indicate 
that the critical behavior is strongly anisotropic, mean- 
ing that correlation lengths diverge with different expo- 
nents in the field direction and perpendicular to the field. 
Specifically, the fluctuations perpendicular to the field 
are gaussian [u± = 1/2) while those parallel to the field 
are not {y\\ = 3/2). Correlations therefore grow faster 
in the parallel direction as T — > T C (E), suggesting an 
analysis of rectangular samples such that the anisotropic 

aspect ratio A = L^^LJ 1 is held fixed [10]. While 
there is some discussion regarding the correct mesoscopic 
model [13], detailed numerical simulations show that the 
exponents cited above are the correct ones [9, 11]. In the 
following we will use only the phenomenological result of 
Leung for the scaling of the order parameter at fixed A: 

S(T i L h L ± )=L/ , ^(tL^,L^L 1 1 ) (6) 

where S refers to 19(1,0). A detailed discussion of the 
subtleties of the FSS analysis for the KLS model and 
precision numerical results can be found in [9]. Fig. 24 
presents our data for the scaled order parameter at 
E = 2; the same data for saturation E can be found 
elsewhere [27]. This data is not intended as a test of the 
mesoscopic model, it merely is meant to indicate that 
the exponents at infinite E are consistent with those at 
finite E, and to determine T c (E = 2.0) = 1.20 (2). The 
data collapse is comparable to that seen at saturation 
E, with the high temperature (upper) branch collapsing 
quite well and the low temperature (lower) branch show- 
ing small, but systematic deviations from scaling. These 
deviations remain unexplained. They are possibly due to 
corrections to scaling or perhaps the asymptotic region 
is only observed very close to T c . 

We have seen how subtle are the issues surrounding the 
transition at / = 0. Perhaps a fruitful way to proceed 
when / =/= is to adopt the technique introduced by 
Caracciolo et al. [28], directly measuring finite volume 



3.0 n 

2.5 

2.0 



1-5 



co 1.0 

0.5 -I 
0.0 



E = 2 



A Xn xV D 



0.0 0.5 1.0 

Ly A (2/3)|T-Tc|/Tc 



1.5 



FIG. 24: Anisotropic scaling plot of S (1, 0) at / = 0. System 
sizes: Squares, 24 x 54; x's, 28 x 86; triangles, 32 x 128. 

correlation lengths for various geometries and volumes. 
In this way we may develop some understanding of how 
to approach the infinite volume limit in a simple way, 
minimizing corrections to scaling which would complicate 
an uninformed analysis. 



IV. CONCLUSIONS AND OUTLOOK 

We have compiled a detailed phase diagram for a sys- 
tem of two species of particles, interacting via attractive 
Ising interactions, and driven into opposite directions by 
an external "electic" field, E. The purpose of our anal- 
ysis was to unify previous studies which were restricted 
either to just one species or to having only excluded vol- 
ume interactions. In the former case, particles order into 
a single strip aligned with the field direction while in the 
second model, the two oppositely driven species form a 
jam in the shape of a transverse strip. We monitor struc- 
ture factors, their fluctuations, and, if necessary, their 
time traces to identify the location and character of the 
transitions. Most of our data are taken at fixed E, vary- 
ing the fraction / of the 'minority' species and temper- 
ature T. At high /, we observe a continuous transition 
from a disordered phase into the transverse strip, as the 
temperature is lowered. Noting a significant system-size 
dependence of the critical line, we invoke a mean-field 
scaling argument [8, 14] which suggests that EL y /T is 
a good scaling variable. This is confirmed quite satis- 
factorily by our data. At smaller / (0.05 < / < 0.10 
for a 40 x 40 system) we observe two transition as T is 
lowered: first into the transverse strip (continuous), and 
then into the parallel strip (first order). And finally at 
the smallest / a single, continuous transition is observed 
into the parallel strip. The junction of the three phases - 
disorder, transverse and parallel strip - appears to be a 
multicritical point. Analyzing data for a range of system 
sizes suggests that its location scales with / /L x , rather 
than /; i.e., the relevant quantity is the number of rows 
which can be formed by the minority particles. Depend- 
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ing on how the thermodynamic limit is approached, the 
multicritical point can shift to / = 0. 

Several projects suggest themselves to extend this 
work. First, an analysis of larger systems at fixed T 
could clear up some questions, especially regarding the 
fate of the 'phase competition' which is observed close to 
the multicritical point. Many of the 'transitions' in the 
40 x 40 system require an analysis of larger lattices in or- 
der to confirm their existence. Second, a look at structure 
factors in the disordered phase is likely to reveal the pres- 
ence of long-range correlations. Since these are known to 
be quite distinct in the KLS [3] and SHZ [17] models, it 
would be interesting to investigate how the crossover oc- 
curs. Such a study would also yield considerable insight 



into the type of noise terms which would have to be added 
to the mean-field equations in order to capture fluctua- 
tions and critical properties accurately. These equations 
would then provide a reliable starting point for an anal- 
ysis of the KLS transition in the presence of a few mi- 
nority charges, in order to understand the true nature of 
the KLS critical point: does it mark the beginning of a 
critical line, or is it a multicritical point? 
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